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ABSTRACT 

Suppose  a  random  vector  X  has  a  multinormal  distribution  with 

k 

covariance  matrix  i  of  the  form  t  =■  I  e  M  ,  where  the  M's  form 

11  1 

1-1 

a  known  complete  orthogonal  set  and  are  the  distinct  unknown 

eigenvalues  of  j:.  The  problem  of  estimation  of  Is  considered 
under  several  plausible  loss  functions.  The  approach  Is  to 
establish  a  duality  relationship:  estimation  of  the  patterned 
covariance  matrix  t  Is  dual  to  simultaneous  estimation  of  scale 
parameters  of  Independent  chi-square  distributions.  This  duality 
allows  simple  estimators  which,  for  example,  Improved  upon  the  MLE 
of  It  also  allows  Improved  estimation  of  tr^.  Examples  are 
given  In  the  case  when  has  equlcorrelated  structure. 

1.  INTRODUCTION  AND  SUMMARY 


Recently  there  has  been  considerable  interest  In  the 
estimation  of  the  covariance  matrix  of  a  multivariate  normal 
distribution.  This  problem  Is  addressed  extensively  In  Stein 
(1975,  1977),  Olkln  and  Selliah  (1977),  Haff  (1977,  1979,  1982), 
and  Dey  and  Srlnivasan  (1985,  1986)  under  plausible  loss  functions. 
However,  there  is  no  work  of  our  kind  available  in  the  literature 
when  the  covariance  matrix  has  an  assumed  structure. 

Suppose  a  random  vector  X  has  a  multinormal  distribution  with 
mean  zero  and  covariance  matrix  t,  which  has  the  form 


E  e  M 

1-1  ‘ ' 


(i.l) 


where  the  9^'s  are  the  distinct  but  unknown  eigenvalues  of  t  and 
the  M^'s  are  a  known  complete  orthogonal  set  of  projection 
matrices.  Such  a  structure  for  t  arises  In  many  practical 
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situations.  A  familiar  example  Is  the  equlcorrelated  case,  that 

Is,  t  “  o  [(1-d)I  +  dJ  ]  where  I  is  the  identity  matrix  and  J 

P  P  P  P 

is  a  p  X  p  matrix  of  I's.  This  is  often  referred  to  as 

intraclass  correlation  structure.  More 

generally,  patterned  covariance  matrices  of  the  form  (1.1)  arise 
naturally  in  variance  component  models.  See  Albert  (1976)  for 
details. 

From  the  classical  viewpoint  one  would  estimate  t  by 
obtaining  the  maximum  likelihood  estimates  of  the  6^'s  using  the 
normality  of  X.  In  this  paper,  however,  we  take  a  decision 
theoretic  approach  for  the  estimation  of  t  using  the  following 
loss  structures: 


and 


L  -  tr(t-t) 

q 

L  (t,t:)  *  trtt  ^  -  log  ill:  ^  |-  p 


(1.2) 

(1.3) 


The  loss  (1.2)  is  the  usual  extension  of  squared  error  loss  (SEL) 
and  the  loss  (1.3)  is  based  on  entropy  measure  of  distance.  Under 

these  losses  the  MLE  is  inadmissible  and  substantial  improvement  is 

available  (see  Table  1).  We  may  show  that  for 

k 

estimators  of  the  form  T  9  m  ,  these  losses  become,  respectively, 

-  ^ 

L(e,e)  -  z  p^Ce  -9  )  , 

i-l  ^ 


L(e,0)  =  ^  "  11  <i-5) 

k 

where  p^  =  rank(M^),  ?  -  ( ,  ...,  0^)  and  E  p^  “  p. 

i-1  ^ 

In  addition,  we  note  the  following  result  essentially  given 
in  Albert  (1  . 
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Theorem  1.1.  Suppose  X  '''  N  (0,^).  Define  Q  =  X'M  X,  where 
-  P  1  i 

k 

1  =  Z  M  ,  M  being  orthogonal  projection  matrices  free  from 

”  1-1  ‘  ‘ 

6  's,  having  rank(M  )  ■  p  ,  i*l,  •••»  k(<  p).  Then  a  necessary  anc 
i  i  1  2  ~ 

sufficient  condition  for  (1)  Q  6  X  •  i”l»  •••>  k,  and 

1  1 


(2)  Q^'s  mutually  independent  Is  j:  =  E 

i“l 

The  equivalence  of  (1.2)  and  (1.4)  and  of  (1.3)  and  (1.5)  for 


estimators  of  the  form  E  6  M  along  with  Theorem  1.1 

i  i  ® 

1»1 

establishes  the  following  duality.  Estimation  of  the  patterned 
covariance  matrix  |  in  (1.1)  under  loss  (1.2)  ((1.3))  is  dual  to 
simultaneous  estimation  of  the  scale  parameters  of  independent 
chi-square  random  variables  under  loss  (1.4)  ((1.5)). 

It  is  to  be  noted  that  in  the  decomposition  of  *  E  6  M  , 

i=l  ^ 

the  0^'s  are  the  distinct  eigenvalues  of  t  with  multiplicity  p^. 
Thus,  for  example,  in  the  equicorrelated  model 


t  =  a^[(l-p)I  +  cj  ],  0  -  o^(l-P),  0  -  0^(1  +  (p-l)P}, 

p  p  1  2 

M,  =  I  “  P  ^  *  p  ^  J  ,  p,  *  p  -  1 ,  and  p  *  1 . 

1  .  p  P  2  pi  2 


In  Section  2,  we  study  the  estimation  of  ^  under  loss  (1.2) 

and  also  tr^  under  SEL.  We  illustrate  for  with  equicorrelated 

2 

structure.  In  this  case,  improved  estimation  of  a  is  discussed. 

2 

Additionally,  Improved  estimation  of  po  is  also  considered. 
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Section  3  is  devoted  to  estimation  of  %  under  loss  (1.3). 
Finally,  in  Section  4,  some  encouraging  numerical  results  are 
given  for  the  equicorrelated  model. 

2.  IMPROVED  ESTIMATION  UNDER  L  LOSS 

-  q  - 

2.1.  Estimation  of 


Here  we  assume  the  conditions  of  Theorem  1.1  and  generalize 

-  k 

•  .  .  .  .  s 

slightly  the  above  discussion  to  estimation  of  ^  =  E  6^  M^  under 

the  loss  ^  ^ 


-  tr(i=  - 


(2.1) 


which,  using  estimators  of  the  form  E  6.  M.,  is  dual  to  the 

i=l  ^ 

s  s  s 

estimation  of  0  =  (6,,..., 6,)  under  the  loss 

1  k 

L(?,e^)  =  E  p.(e®  -  e!)^. 

i=i^  ^  ^ 


(2.2) 


For  the  estimation  of  s  =  1;  for  the  estimation  of  the  precision 

“1  2 
matrix  ^  ,  s  =  -1.  Let  Q  =  (Q-,...,Q,  )  where  Q.  'v  e.x  and  are 

X.  iC  X  X  o  • 

0  s  ’-, 

independent.  Suppose  6  (Q)  is  an  estimator  of  9  given  componentwise 

Os  ^  s 

as  6£(Q)  =  ^  “  l,...,k.  For  example,  a^  =  p^  gives  the  MLE  of 


Now  define 


(1)  r.  ^=E(QT|e.  =  l)/E(Q?|e.  =  1) 

i,a,6  i'  1  i'  1 

=2“"®  r((p.+2a)/2)/r((p.+26)/2). 

k  a . 

(2)  For  a  =  (ai,...,a,  ),  v  *  H  E(Q.’’|9.  =  1) 

1  k  a  1  '  1 


k  a . 

=  n  2  ^r((p.+2a^)/2)/r(p^/2). 


Note  that  a.  =  r.  „  gives  the  best  invariant  estimator  of  6.  under 

XXjSj^S  X 

squared  error  loss. 


Thus  If.  =  Lp.^Q®  is  the  maximum  likelihood  estimator  of  t  and 

Is  -  s  .  ,  .  .  . 

4-  =  Zr.  ^  Q.  IS  the  estimator  obtained  by  combining  the  best 
vj  1)S)^S  1 

•  •  »  s  .  is 

invariant  estimators  of  6^.  The  following  lemma  shows  that 

dominates  if,  under  the  risk  criterion- 

Lemma  2.1.  1  ^4^4^)  ^ 

Proof.  Immediate  from  the  duality  (2.1),  (2.2)  and  the  fact  that 
o  9=Q-  dominates  p.  Q.  in  estimating  0.. 

In  view  of  Lemma  2.1,  it  is  sufficient  to  find  estimators  which 
improve  |q.  We  have  the  following  theorem  whose  proof  is  a  special 
case  of  a  result  in  Dey  and  Gelfand  (1987). 

Theorem  2.1.  Consider  the  estimator  6(Q)  =  (6^ (Q) , . . . , 6  (Q) ) 
given  componentwise  as 


6.(Q)=6?(Q)  =b(  n  i  =  1,...  ,k(>  2), 

^  j=l  J 


(2.4) 


where  6.(Q)  =  r.  ,  Q..  Then  provided  all  expectations  exist,  6(Q) 

1  l^Sj^Sl 

dominates  6*^(Q)  under  loss  (2.2)  if 

n  .  K  ^  9  .4^2) 

where  d^^^  =  min(p.|d.j)  and  d^^^  =  max(p. )  with 

r  1  I  1  I  r 


‘^i  ’^i  ,s  ,2s^i  ,  (k+l)s/k,s/k 


-  1,  i  =  1, . . . ,k. 


In  view  of  Theorem  2.1,  it  follows  that  under  the  loss  (2.1),  an 
improved  estimator  of  is  given  as 


i®  =  Z  6.(Q)M.  =  |!  +  b(  n  Q.)®^*'  I  . 
i=l  ^  ^  0  j=l  J  P 


(2.5) 


Remark  2.1.  Theorem  2.1  requires  only  that  follow  a  distribu¬ 
tion  with  6^  as  scale  parameter.  In  this  setting  Dey  and  Gelfand 
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(1987)  offer  more  general  results  than  Theorem  2.1  including  estimators 

which  provide  maximum  improvement  along  a  ray  determined  by  a  specified 

vector  0Q.  These  estimators  shrink  (expand)  differently  on  each 

coordinate.  Using  the  fact  that  the  are  distributed  as  multiples 

of  chi-square  random  variables  enables  Klonecki  and  Zontek  (1987)  to 

provide  necessary  and  sufficient  conditions  for  the  existence  of  an 

estimator  of  the  form  (2.4)  to  dominate  6^  given  componentwise  as 
0  s 

6^  =  for  any  specified  a^,  i  =  l,...,k. 

Example :  Suppose  \  has  equicorrelated  structure.  In  this  case 
2  2 

9^  =  0  (1  -  p),  0^  =  o  {1  +(p  -  l)p},  are  the  distinct  eigenvalues  of 
The  best  invariant  estimate  of  6^  is  6?(Q)  =  (p^  +2)  ^  with 
p^  =  p  -  1  and  p^  =  1,  a  class  of  improved  estimators  of  9  =  (0^,02) 
is  given  componentwise  as 


—  1  ^  k 

6.(Q)  =  (p.  +  2)  "q  +  b(  n  Q  y,  i  =  1,2, 

^  j  =  l  J 


(2.6) 


where  0  <  b  <  2v,  d^  ^  Vv,  d^^^  and  the  upper  bound  on  b  simplifies  to 
?  1 

2r(p/2)/3(p-i)r((p+i)/2)A’. 

The  corresponding  improved  estimator  of  \  has  the  simple  form 
f  =  Jq  +  b(Q^Q2)^  I  .  (2.7) 


Remark  2.2.  Improved  estimation  of  0  ^  =  ( 0^^ , .  .  .  ,0^^^)  ,  hence  of 

^  follows  directly  from  Theorem  2.1.  We  only  need  the  existence 

of  appropriate  reciprocal  moments  of  the  Q^.  Unfortunately,  in  the 

—a 

equicorrelated  case  EQ^  does  not  exist  for  a  ^  1  and,  hence,  our 
approach  does  not  provide  a  dominating  estimator. 
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2.2.  Estimation  of  trace  of  !! 

Consider  now  estimation  of  the  trace  of  |  under  the  SEL  given 


L(a,tr|i)  =  (a  -  tr^)  . 


(2.8) 


Since  tri  =  Z  p.0.  our  duality  converts  estimation  of  the  trace  to 
i=l  ^ 

estimation  of  a  linear  combination  of  chi-square  scale  parameters. 

The  following  theorem  gives  a  class  of  admissible  estimators  of 
k 

Z  £.6.  if  £.  >  0. 
i=l  "  " 


Theorem  2.2.  If  ct.  >  0,  i  =  l,...,k,  known  then  under  the  loss 


(2.8) 


Z  £.a.  ^ 

V“i 

1=1 


(2.9) 


is  admissible  for  Z  £.0.. 

i=l  "  " 

Proof.  Consider  the  subset  of  the  parameter  space 

C  =  {(0^,.  ..  ,0j^)  :  0^  =  a^0,  >  0,  i  =  l,...,k}. 

2  ^ 

Then  Q./o.  ©x  .  i  =  l,...,k.  Thus,  on  C  Z  Q./a.  is  sufficient  for 

11  P .  .,11 

k  1  2 

0  and  Z  Q./a.  0x  Thus,  by  a  theorem  of  Karlin  (1958), 
k  i=l  ^  ^  p' 

Z  Q./a. 

■=1  ^  ’■ 

— - s -  is  admissible  for  0  under  SEL.  Then  on  C,  6  (Q) 

p+2  a 


=  {(  Z  £.a.)  Z  Q./a.}/(p+2)  is  admissible  for  Z  £.a.6  =  Z  £.0.. 


1  1  .  ,  1  1 
1=1  1=1 


Suppose  '5^(Q)  is  inadmissible  for  Z  £.0.,  Then  there  exists  '^^(Q) 

i=l 

which  dominates  5^(Q).  But  '^^(Q)  being  admissible  on  C  implies 

5*(Q)  =  6  (Q)  a.s. 
a  a 
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Note  that  since  >  0  in  the  expression  for  tr^  Theorem  2.2  applies. 

Remark  2.3.  Consider  the  equicorre lated  structure.  As  a  special 

case,  =  1  gives  X'X  admissible  for  trace  |  and  hence  X  X/(p+2) 

2 

admissible  for  o  .  Similarly  a.  =  p_/{p.+2)  implies 
2  1-  i  1 

Z  p. /(p.+2) 


i  =  l 


p+2 


Z  p.Q./(p.+2)  is  admissible  for  o  ,  i.e.,  an  appropriate  linear 
i=l  ^  ^  ^ 


combination  of  the  componentwise  best  invariant  estimator  is  admissible. 

Now  we  will  demonstrate  a  general  method  for  improving  on  a  linear 
estimator  of  a  linear  combination.  The  improved  estimators  are  nonlinear, 
and  may  shrink  or  expand  the  given  linear  estimator.  Work  of  Das  Gupta 
(1986),  Dey  and  Gelfand  (1987),  and  Klonecki  and  Zontek  (1987)  is  relevant 
here.  A  general  result  is: 

Theorem  2.3.  Provided  expectations  exist,  an  estimator  of  tr^ 
k  k 

of  the  form  6.  =  Z  Jl.Q.  (l.  =  1  yields  Z  Q.  =  X’X  the  MLE ,  which  is 
0  1^1  1  ^  .  ,1 
i.=l  1=1 

also  UMVUE)  is  dominated  by 

k  r . 


6  =  6  +  c  n  Q . 

0  j=l  J 


J 


(2.10) 


where  r.>^0,  Zr.  =  l,  if  and  only  if  either 
J  j  =  l  J 

(i)  ^  0 ,  r^  =  0  if  d^  =  0  an  1  c  >  0  sufficiently  small, 

(ii)  .1  =  0  if  d^  =  0  and  c  <  0  sufficiently  large, 

where  d.  =  (l-£.)p./2  -  2..r.,  i  =  l,...,k  and  d,,  =  min  d.  ,  d.,  ,, 

1  11  11  (1)  1  (k) 

=  max  d . . 

1 

Proof.  The  risk  difference  between  (2.10)  and  6^  is 


k  2r. 


k  r. 


A(0)  =  c  V.,  n  e.  ^  -  cv  n  e.^  z  e.d.  . 


2r  j 


j=i  J  j  =  l  J  J 
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Suppose  (i)  holds  (the  proof  for  (ii)  is  similar)  and  that  all 


i  such  that  r.  =  d.  =  0  have  been  deleted  in  A(0).  Then 
J  J 

k 

E  e.d. 

\  i=l  J  J 

^2r  k  r . 


k  2r . 

A(0)  =  cv.  n  6.  ^[c  - 
2^  i=l  J 


]. 


n  0. 

j=i  J 


J 


r .  r . 

But  E0.d./n6.'^  >  n(d./r.)  .  Hence,  A(0)  <  0,  V0 ,  if  c  is  positive 

J  J  J  -  J  J 

and  sufficiently  small. 

Remark  2.4.  In  particular,  when  =  1  any  set  of  nonnegative 

r.  such  that  Er.  =  1  and  at  le  ist  two  r.  differ  from  zero  will  work. 
11  1 

Here  c  <  0  so  that  the  dominating  estimator  is  a  shrinker.  Since 

tri  >  0,  6  =  max(6  ,0)  will  dominate  6  (using  a  lemma  of  Klotz, 

r,c  r,c  ^ 

Milton  and  Zacks  (1969),  p.  1394). 

Example .  Again  consider  the  equicorrelated  structure.  Clearly 
2 

trj  =  po  .  Thus,  using  (2.10),  we  can  explicitly  dominate  X’X  in  esti- 
2  .  2 

mating  pa  ,  hence  X'X/p  in  estimating  o  by  nonlinear  estimators.  For 


instance,  the  estimator 

5*  =  X'X/p  +  6(0^0^)^ 

2  2 

dominates  the  MLE  under  loss  function  (S  -  a  )  if 


(2.11) 


-  4r(p/2)/pr(£^)  A  <  b  <  0. 

Remark  2.5.  If  we  attempt  to  apply  Theorem  2.3  to  linear 

estimators  of  the  form  (2.9),  we  will  discover  that  all  d^'s  are 

equal  to  zero.  In  particular,  in  the  above  example,  we  cannot 

dominate  X  X/(p+2)  in  estimating  a  . 

Continuing  with  our  example,  suppose  p  ^  0  and  consider  estimation 
2 

of  pa  which  may  be  viewed  as  variance  component.  (See  Gelfand  and  Dey 
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(1988)  for  more  detailed  discussion  of  improved  estimation  of  variance 


components.)  Since  po  =  (62~5^)/p,  we  consider  the  estimator 
5^  =  p  ^(^2^2  ~  ^1*^1^ 

where  a^  =  (p^  +  2c  '  \  i  =  1,2,  with  0  ^  £  1-  For  example,  taking 

=  0,  i  -  1,2,  5^  becomes  the  MLE  of  oo^  and  taking  =  1,  6^  is  formed 


from 


the  best  invariant  estimator  of  6^,  under  SEL,  i  =  1,2.  A  class 

2  . 


of  improved  estimators  of  pa  ,  is  given  as 


6=6°+  b(a^Q^)"(a2Q2)^  (2.12) 

using  Theorem  2.3,  provided  r  >  0  can  be  chosen  such  that  either  (i) 
r  <  min(ej^,  1  -  e^)  whence  b  must  be  positive  and  sufficiently  small 
or  (ii)  r  >  max(£j^,  1  -  whence  -b  must  be  positive  and  small.  In 
fact,  we  would  use  6^- 

Remark  2.6.  Note  that  Theorem  2.2  does  not  provide  admissible 
2 

estimators  of  pr  since  1.  <  0. 

1 

3.  IMPROVED  ESTIMATION  UNDER  L  LOSS 

-  c  - 

In  this  section  we  consider  the  estimation  of  patterned  ^  under 

loss  (1.3).  Using  the  aforementioned  duality  for  estimators  of  the 
k  . 

form  Z  9.M.,  we  convert  this  problem  to  simultaneous  estimation  of 
i=l  ^  ^ 

the  eigenvalues  6  =  (9^,...,6j^)  of  ^  under  the  loss  (1.5).  In  fact, 

s  s 

s 

our  results  can  be  extended  to  the  estimation  of  9  =  (6,  . 6,  ) 

1  k 

and  hence  the  estimation  of  .  However,  it  is  not  clear  how  to  apply 
loss  (1.3)  in  estimating  tr|  which  is  not  a  scale  parameter. 

Our  approach  is  that  of  Berger  (1980)  and  Dey,  Ghosh  and  Srinivasan 
(1987).  Assume  are  independent  Gamna(a^,n^)  (a^  >  0,  n.  >  0) 
random  variables,  having  density 
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ot.  a._,  -Q.n. 

f(Q.|n-)  =  n-  Q.^  e  ^  ^/r(a.) 

i'  111  1 


(3.1) 


In  our  case,  we  have  a.  =  p./2,  n-  =  (26.)  ,  whence  the  loss  (1.5) 

111  1 

corresponds  to  simultaneous  estimation  of  the  that  is, 


-1 

L(6,n  )  =  2  p.[6.n.  -  log(6.ri.)  -  l]. 


Since  the  MLE  of  9.  is  6.(Q)  =  Q./p.,  the  MLE  and  the  unbiased  estimator 
1  1  11 

of  6.  =  r\.^l2  is  6?(Q)  =  Q./2a.,  i  =  l,...,k.  From  Dey ,  Ghosh  and 
11  1  11 

Srinivasan  (1987),  it  follows  that  6*^(Q)  =  (6^(Q)  , .  .  .  ,  6^(Q) )  is  admissible 
for  k  =  2  if  min(aj^,aT)  >  4.  To  seek  a  dominating  estimator,  we  require 
k  >  2. 

Now  assuming  the  conditions  in  Lemma  1  of  Berger  (1980),  it 
follows  that  if 


6(Q)  =  5  (Q)  +  i(Q) 


(3.2) 


is  a  rival  estimator,  the  risk  difference  is 
A(e)  =  R(6,e)  -  R(6°,6)  =  EgAplQ), 

where  A^(Q)  is  the  unbiased  estimate  of  the  risk  difference  given  as 


Aq(Q)  =  I  p.[^J^^^(Q)  +  (a.-DA.  (Q)/Q^  -  log(  1+a  .  i  .  (Q) /Q.  )  ] 
i=l 


with  (Q)  =  3(f!^(Q)/3Q^.  Defining  1;^(Q)  =  Q^4.^(Q)  ,  one  gets 


A  (Q)  =  Z  p.[Q.i';.  (Q)  +  a.ii;.(Q)  -  log(  1  +  a.i^.(Q))]. 

0  .  ,111  11  11 
1  =  1 


(3.  3) 


In  order  to  obtain  an  improved  estimator  6(Q),  it  is  sufficient 


to  find  a  solution  ^^(Q)  ^  0  with  strict  inequality  for  some  set  of  Q. 
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The  following  Lheorem  gives  a  class  of  dominating  shrinkage 
estimators . 

2 

Theorem  3.1.  Suppose  S  =  Z  log  (Qj^/2).  Consider  ar  estimator 
5(Q)  =  (5  j^(Q)  , .  . .  ,  6^(Q)  )  given  componentwise  as 

6^(Q)  =  pT^Q.  -  Q^T(S)log(Q./2)/2{b+S),  i  =  1 . k,  (3.4) 

2  *2 

with  b  >  (5.76)(k-2)  /p  ,  where  p*  =  max  p^  and  t(S)  is  a  function 
satis  fying 

(i)  0  <  t(S)  <  4.8(k-2)/p*^ 

(ii)  t(S)  f  in  S  and 

(iii)  E[t'(S)]  <  CO.  (3.5) 

Then  6(Q)  dominates  6*^(Q)  for  k  ^  3,  in  terms  of  risk. 

Proof.  The  argument  is  similar  to  that  of  Theorem  3.1  of  Dey , 

Ghosh  and  Srinivasan  (1987). 

Remark  3.1.  Using  Theorem  3.2  and  3.3  of  Dey,  Ghosh  and  Srinivasan 
(1987),  adaptive  estimators  and  trimmed  shrinkage  estimators  of  6  can 
be  obtained  as  well. 

In  concluding  this  section,  we  observe  another  illustration  of  our 
duality  relationship.  Improved  estimation  of  patterned  |  under  the 
scale  invariant  loss 

L  (|,Z)  =  tr(^^"^  -  I)^  (3.6) 

^  k  „ 

using  estimators  of  the  form  Z  9.M.,  converts  to  simultaneous 

i=l  ^  ^ 

estimation  of  S  =  (e^,...,6j^)  under  loss 
k  . 

L(e,e)  =  z  p.(9./e.  -i)  . 

t  1  1 


(3.7) 
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Again,  Berger's  approach  yields  a  differential  inequality  (the 
only  difference  will  be  the  presence  of  the  weights,  p^)  whose 
solution  leads  to  dominating  estimators  similar  to  those  in  (3.4). 
Details  are  omitted. 

4.  NUMERICAL  RESULTS 

2 

To  study  the  performance  of  the  MLE  =  E  p.  Q. , 

M  .,11 

2 

^  ^  —  1  ..K  L 

Iq  =  I  (p^  +  2)  and  ^  ^  +  b(Qj^Q2)^  1  ,  we  calculate  risks 

i=l  ^  P 

for  different  values  of  p  and  p  in  the  equicorrelated  structure.  We 

2 

took  o  =1  and  b  =  T  (p/2 ) /3(p-l)r(  (p+1 ) /2) /ir ,  which  is  the  midpoint 
of  the  allowable  range.  We  then  computed  the  percentage  improvements 
for  selected  values  of  p  and  p  (-(p-1)  ^  ^  P  ^  The  improvements  of 

I  over  the  MLE  are  substantial.  While  the  percentage  improvements 
in  risk  of  |  over  are  small,  the  simplicity  of  \  encourages  its 


use. 
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TABLE  1 


PERCENTAGE  IMPROVEMENTS  OVER  AND  Z„ 

M  0 


p 

R(L)  -  R(j) 

PIl  =  - -  X  100 

R(i„) 

R(L)  -  R(i) 

PI2  =  - -  X  100 

’"to' 

-.75 

67.94 

3.83 

-.50 

68.43 

5.31 

-.25 

68.78 

6.  35 

0 

68.92 

6.75 

.25 

68.78 

6.35 

.50 

68.44 

5.31 

.75 

67.94 

3.83 

P  =  6 

0 

49.25 

3.11 

.25 

63.77 

2.45 

.50 

66.42 

1.50 

.75 

66.85 

.85 

p  *  10 

0 

43.50 

1.87 

.25 

64.66 

1.23 

.50 

66.50 

.  70 

.  75 

66.75 

.  39 
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